Collective Modes and Skyrmion Excitations in Graphene SU(4) Quantum Hall 

Ferromagnets 



o 
o 

(N 

wo: 
< 



Kun Yang 

National High Magnetic Field Laboratory and Department of Physics, 
Florida State University, Tallahassee, FL 32306, USA 

S. Das Sarma 
Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, MA 20742, USA 



A. H. MacDonald 

Department of Physics, University of Texas, Austin, TX 18112, 
(Dated: February 6, 2008) 



USA 



Graphene exhibits quantum Hall ferromagnetism in which an approximate 517(4) symmetry in- 
volving spin and valley degrees of freedom is spontaneously broken. We construct a set of integer 
and fractional quantum Hall states that break the 517(4) spin/valley symmetry, and study their 
neutral and charged excitations. Several properties of these ferromagnets can be evaluated analyt- 
ically in the SU(4) symmetric limit, including the full collective mode spectrum at integer fillings. 
By constructing explicit wave functions we show that the lowest energy skyrmion states carry charge 
±1 for any integer filling, and that skyrmions are the lowest energy charged excitations for graphene 
Landau level index |n| < 3. We also show that the skyrmion lattice states which occur near integer 
filling factors support four gapless collective mode branches in the presence of full 5(7(4) symmetry. 
Comparisons are made with the more familiar SU (2) quantum Hall ferromagnets studied previously. 
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I. INTRODUCTION 

Recent experimental worki has established graphene 
as a new two-dimensional (2D) electron system with 
linear Dirac-like energy-band dispersion. An impor- 
tant aspect of graphene physics is the two-fold valley 
degeneracy of its low-energy band structure which, in 
combination with the usual spin doublet, implies band 
eigenstates with degeneracy N = 4. In a strong mag- 
netic field, the four-fold degeneracy of Landau levels in 
graphene has been clearly observed in recent quantum 
Hall measurements This property is expected^ to 
lead to an intriguing interplay between interaction and 
quantum Hall physics by introducing a rich variety of 
new states with different types of spontaneous symmetry 
breaking, and new varieties of low- lying collective modes. 
In this paper we explore some of the theoretical possibili- 
ties for quantum Hall ground states and collective modes 
that follow from the enlarged (TV = 4) Landau level de- 
generacy of graphene, emphasizing similarities and dif- 
ferences compared to the well-studied N = 2, spin-only 
circumstance^ relevant to GaAs 2D systems. 

One of the most theoretically interesting 8 - 9,10 and phe- 
nomenologically significant*^ consequences of quantum 
Hall ferromagnetism is the presence of a finite density of 
skyrmions in the ground state near integer filling factors. 
A dense system of skyrmions constitutes an emergent 
set of low-energy degrees of freedom that qualitatively 
altersii NMR, optical, thermal and transport properties 
of the two-dimensional electron system. In the N = 2 
case, dense skyrmion systems occur only near Landau 
level filling factor v = 1. We predict that in graphene 
dense skyrmion states&Siii occur near many integer fill- 



ing factors, and that they have a larger number of gap- 
less collective modes than for N = For Coulombic 
electron-electron interactions, skyrmion lattice states oc- 
cur onlyi^ near v = 1 in the £{7(2) case. For graphene 
we predict skyrmion lattice states near v = ± I for all 
/ < 14 except for I = 2, 6, 10 when all Landau-level mul- 
tiplets are either full or empty, and that skyrmion lattice 
states in graphene have four branches of gapless collec- 
tive modes at all these filling factors in the absence of 
symmetry-breaking perturbations . 

The N — 4 internal degrees-of-freedom present in 
graphene take on particular significance because of the 
very strong Coulomb interaction energy scale (estimated 
to be more than 1000K at A5T^ several times larger than 
in GaAs) which will help make all the physics that follows 
from interaction-driven, spontaneous symmetry breaking 
(i.e. phases which break the internal symmetry associ- 
ated with the degeneracy N) more accessible experimen- 
tally. The point of departure for our analysis will be an 
SU(N — 4) symmetric Hamiltonian which allows us to 
obtain a number of exact results. This highly symmet- 
ric model is believed to be a good approximation of the 
full Hamiltonian of graphene For relatively low car- 
rier densities the largest symmetry breaking term in the 
Hamiltonian will be Zeeman coupling which favors spin 
alignment along the field direction. We will discuss some 
of the important consequences of this term. 

We remark that the possibility of a spontaneously bro- 
ken valley symmetry is akin to orbital ordering, with 
the novel twist that the orbital ordering is in momentum 
space, rather than in real space. The anticipated exotic 
orbital order in graphene should co-exist with spin order 
making graphene (along with systems such as superfluid 
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He-3, manganites, and related systems where real space 
orbital order apparently coexists with spin order) an in- 
teresting system to study novel quantum phases with in- 
terplay between spin and orbital order. The momentum 
space orbital ordering discussed here may also imply ob- 
servable Jahn- Teller coupling effects in graphene, but we 
do not explore this idea further in this paper. 

The consequences of Si/ (4) or more generally, SU(N) 
symmetry in quantum Hall physics have been discussed 
previously by Arovas and co-workers^ in the context of 
Silicon systems which also have additional (approximate) 
valley degeneracies, and by Ezawa and co-workersi&ii in 
the context of bilayer quantum Hall systems where the 
additional degeneracy comes from the which layer degree- 
of-freedom. Although both groups used non-linear sigma 
model descriptions they reached different conclusions on 
the properties of skyrmion excitations in SU (N) quan- 
tum Hall ferromagnets, and in particular on the mini- 
mum charge a skyrmion can carry. By explicitly con- 
structing the microscopic wave functions for collective 
modes and skyrmions, we demonstrate that the lowest- 
energy skyrmion excitations of the system carry charge 
±1. This analysis also allows us to enumerate the in- 
ternal degrees of freedom associated with an individual 
skrymion and predict the number of gapless collective 
modes present in SU(A) skyrmion lattice states. 

Our paper is organized as follows. In Section II we 
identify exact eigenstates of SU(N) quantum Hall ferro- 
magnets at integer filling factors which have broken sym- 
metry. We argue that these states are the ground states 
for any physically sensible repulsive interaction between 
the particles, and obtain exact results for their elemen- 
tary particle-hole excitations. We find that in the absence 
of symmetry breaking perturbations there are M (N—M) 
gapless collective modes when M of the N members of 
a Landau-level multiplet are occupied. In Section III 
we discuss the properties of skyrmions in graphene. We 
find that for Coulombic electron-electron interactions, 
skyrmions are much more robust in systems like graphene 
with a Dirac band structure, than for the more familiar 
parabolic band systems. We predict that dense-skyrmion 
systems will be ubiquitous in the quantum Hall regime 
of graphene, and that at even integer fillings skyrmion 
lattice states will have four branches of gapless collective 
modes, two with k 1 dispersion, one with fc 3 / 2 dispersion 
and one with k 2 dispersion. Finally in Section IV we 
briefly present some elementary considerations on broken 
symmetry states in graphene at fractional filling factors. 
We conclude in Section V with a brief summary of our 
results and some final comments. 



II. SYMMETRY BREAKING GROUND STATES 
AND COLLECTIVE MODE SPECTRA 

For simplicity we ignore disorder and mixing between 
different Landau levels. To keep the discussion general, 
we assume electrons have N internal states and that the 



electron-electron interaction is independent of these in- 
ternal states, so that the Hamiltonian is SU (N) symmet- 
ric. Thus the total degeneracy of each Landau level (LL), 
including the internal degeneracy, is N x 7V$, where 7V$ 
is the number of flux quanta enclosed in the system. We 
start by considering the case where the filling factor of 
the valence LL (with LL index n) is an integer: 

u n = N e /N§ =M<N, (1) 

where N e is the number of electrons occupying the va- 
lence LL. We note that in the case of graphene N = 4 
and the four-fold degenerate n = LL is half-filled (cor- 
responding to Vq = 2) at zero doping; thus the Hall con- 
ductance is 

a xy = ve 2 /h (2) 

with 

v = 4n-2 + i/„. (3) 

In the absence of interactions the quantum Hall effect 
occurs only when each of the four-fold degenerate Landau 
levels is completely full (y n — N — 4) or completely 
empty (v n = 0) and the quantized Hall conductance is 

o~xy = (in + 2)e 2 /h, (4) 

where n is the index of the highest filled Landau level. 
Quantum Hall Ferromagnetism will£ lead to quantum 
Hall pleateaus characterized by intermediate integers and 
by fractions. In the following we neglect interaction in- 
duced mixing between orbitals with different n. For dif- 
ferent values of n the properties of the SU (N) quantum 
Hall ferromagnet are different. In the following the de- 
pendence on n is sometimes left implict to avoid nota- 
tional clutter. 

Because of the SU(N) symmetry of the system, the fol- 
lowing single Slater determinant state (in which only elec- 
tronic states in the valence LL are explicitly described) 
is an exact eigen state of the Hamiltonian: 

i*o)= n ii c U o >- ( 5 ) 

l<a<M k 

Here c 1 is the electron creation operator, |0) is the vac- 
uum state, <7 is the index of the internal state that runs 
from 1 to N, k is an intra LL orbital index; for exam- 
ple in the Landau gauge it is the wave vector along the 
plane wave direction, while in the symmetric gauge it is 
the angular momentum quantum number. 

For a broad class of repulsive interactions, we expect 
|\&o) to be the exact ground state of the system; for the 
case of v n = M = 1 this can be proved rigorously for 
short-range repulsion (see below). In this paper we use 
this assumption as the starting point of our discussion. 
Obviously, the ground state |^o) breaks the SU(N) sym- 
metry spontaneously since an SU (N) rotation transforms 
the M spontaneously chosen occupied single electron or- 
bitals to another different set. It therefore represents the 
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ground state of an SU(N) ferromagnet, and is expected 
to support ferromagnetic spin waves. In the following we 
construct the exact single spin wave states and determine 
their spectrum in a manner similar to that of Kallin and 
Halperin>i& Consider the following Landau gauge states, 
with y in the direction of the plane waves: 




— I kx i ky c 



C fe' + fe y ,<T 2 C fc'^ll^0)- 



(6) 



These spin-wave states are labeled by a two-dimensional 
wavevector k and two internal indices: 1 < <j\ < M and 
M < (72 < N. £ — y^hc/eB is the magnetic length. 
Physically they can be understood as single particle-hole 
states formed from the <T2 and o~\ internal states, i.e. as 
magnetic excitons. Since there are M choices for u\ and 
N — M choices for <72, the total number of these spin wave 
modes is M(N — M). It follows from translational in- 
variance that k is a good quantum number, from SU (N) 
invariance that excitons with distinct (02, ai) labels are 
uncoupled, and that |k) CTljCr2 is therefore an exact eigen- 
state of the Hamiltonian. The exact SU(N) quantum 
Hall ferromagnet magnetoexciton dispersion relation is 

E(k) = (k|V|k) - (* |^|*o> 
1 
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qV(q)[F n (q)] 2 e-* ' ' 2 [l - J (qke 2 )]dq,(7) 



where V(q) is the Fourier transform of electron-electron 
interaction, Jo is the Bessel function, and 

F n (q) = i[L|„|( g 2 £ 2 /2) + L|„|_ 1 (^ 2 /2)] (8) 

is the Landau level structure factor appropriate for Dirac 
fermions. 5 ^ (In Eq. (jHJ) L n (x) is a Laguere polynomial. 
The dependence of graphene quantum Hall ferromagnet 
properties on n enters only through F n .) The second 
term in Eq. is due to particle-hole attraction^ and 
vanishes for k — > 00 where E{k) approaches the energy 
of a widely separated electron-hole pair. For the case of 
a 1/r Coulomb interaction and n = (Fq(x) = 1), we 
have 



E(k) 



e 



1 _ e - k2e2 ^I Q (k 2 £ 2 /4) 



(9) 



where e is the effective dielectric constant and Iq is the 
modified Bessel function. In the long-wave length limit 
E(k) oc k 2 as expected for ferromagnetic spin waves. 

It is worth pointing out here that among the various 
sources of perturbations that break the SU(A) symmetry 
in graphene, the simplest but most important one is the 
Zeeman splitting of electron spin: 



(10) 



where g is the electron spin g factor and [13 is the Bohr 
magneton. Due to the fact that the total spin along the 



magnetic field direction Sf ot commutes with the SU(4) 
invariant Hamiltonian, it only shifts the energies of in- 
dividual eigenstates without changing the states them- 
selves. In the presence of Hz electrons will choose to 
occupy spin-up states in the ground state Eq. (J5J until 
all such states are filled; the internal states as labeled by 
cr are thus eigen states of S z . The single spin- wave states 
take the same form as in Eq. @; the spectrum remains 
the same as Eq. (J7J if u\ and 02 are of the same spin 
orientation, while it gets shifted to E(k) + gusB if they 
are of opposite spin orientations. 

Returning to the SU (N) symmetric case, wc note that 
the nature of the symmetry breaking as realized in the 
ground state (J5J, as well as the resultant gapless spin 
wave excitations, may be understood from the following 
group-theoretical analysis. The state JSJ, while not in- 
variant under a general SU (N) transformation, is invari- 
ant under a subgroup of SU(N), SU(M) x SU(N - M). 
Physically the SU(M) x SU(N-M) subgroup represents 
SU(M) transformations among the M occupied levels, 
and SU(N — M) transformations among the N — M un- 
occupied levels; these transformations do not change the 
state ©. Thus the order parameter as represented by 
the symmetry breaking state (JSJ forms a coset space 
of SU(N)/SU(M) x SU(N - M) x 17(1) (where the 
last U(l) represents an overall phase difference between 
the occupied and unoccupied levels), or equivalently, 
U(N)/U(M) x U(N - M). Since a U(N) transfor- 
mation is parameterized by N 2 parameters (its num- 
ber of generators), we find the dimensionality (or the 
number of independent fields) of the coset space to be 
iV 2 - M 2 - (N - A/) 2 = 2M(N - M). Because we are 
dealing with a ferromagnetic state (whose effective action 
contains a dynamic term with a single time derivative, see 
below), half of these fields are the conjugate momenta of 
the other half; we thus expect M(N — M) independent 
spin-wave modes, in agreement with the microscopic con- 
struction®. Very similar analyses were performed in 
Refs. Il5llla that lead to the same conclusion. 

The analysis above suggests the following matrix 
parametrization of the ferromagnetic order parameter, 
appropriate for the symmetry breaking corresponding to 
U(N)/U(M) x U(N - M)±&: 

Q(r,t) = U\T,t)SU(r,t), (11) 

where U (r, t) is the (space-time dependent) N x N U(N) 
transformation matrix, Sij oc &y is a diagonal matrix 
with S u = 1 for < i < M and S u = -1 for M <i< N. 
The N x N matrix field Q(r,t) is the order parameter. 
Ezawa and co-workers^ use a different, but presumably 
equivalent, parametrization of the order parameter. 

In terms of the matrix field Q(r,t), the long- wave 
length, low-energy effective action of the system takes 
the formic 

S[Q{r,t)\ = J dtdv {atiA(Q) ■ d t Q 

+ 7P.tr(VQ) • (VQ) + •••}, (12) 
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where the first term is the Berry phase termi^ that en- 
codes the commutation relations between different com- 
ponents of the order parameter field, and the second term 
describes the energy cost when the order parameter is 
non-uniform. The choice of the prefactor 1/4 (instead 
of 1/2) for the second term is for later convenience as 
it compensates for the fact that the trace of the square 
of Pauli matrices is two. Terms involving higher orders 
of either time or spatial derivatives, as well as possible 
symmetry breaking perturbations are represented by • • • . 
Thanks to the knowledge of the exact ground state, the 
order parameter stiffness p s may be determined exactly 
in a manner similar to that of Ref. 0: We first con- 
struct a state by performing a slow, r-dependent SU (N) 
rotation on l^o), then project it to the appropriate LL, 
calculate its energy, and perform a gradient expansion of 
the energy functional. The result is 
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2 p -q 2 i 2 /2 



dq. 



(13) 



and for the case of n = and 1/r Coulomb interaction, 
p s = e 2 /(16^/2ne£)}^ p s is independent of both N and 
M; Eq. JTjJJ is identical to the N = 2 and v n = M = 1 
p s expression derived in earlier work 8-19 . This finding 
is not surprising since any infinitesimal rotation in the 
U(N)/U(M)xU{N-M) subgroup of SU(N) can be de- 
composed into combinations of SU (2) rotations between 
occupied and unoccupied levels. It is easy to verify that 
the action l|12f) reproduces the spin wave spectrum in 
the long wave-length limit. 



III. SINGLE SKYRMION STATES AND 
COLLECTIVE MODES OF SKYRMION 
LATTICES 

The matrix field Q supports topologically non-trivial 
spatial configurations, which can be parameterized by 
an integer- valued topological quantum number called the 
Pontryagon index: 

q = [ d 2 ve^tx[Qd^Qd v Q], (14) 

107T J 

where e Miy is the antisymmetric tensor. Field configura- 
tions with nonzero q carry topological charge and are 
called skyrmions. As in the N = 2, v n = M = 1 
case^SiiS quantum Hall ferromagnet skyrmions also carry 
an electric charge that is equal to its topological charge. 
It follows from the above considerations, as concluded 
in earlier work^ that skyrmions with charge ±1 exist 
within the lowest Landau level. This remarkable prop- 
erty implies that when skyrmions are the lowest energy 



charged excitations, they will be present* 



8,9,10,20 



the 



ground state of the system when v is close to, but not 
equal to, an integer. Skyrmions thus appear as emergent 
low-energy degrees of freedom and influence all observ- 
able properties. 
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1/2 
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1/2 
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11/16 


7/32 


3/4 


7/8 




(0.6875) 


(0.2188) 


(0.75) 


(0.875) 


2 


145/256 


169/512 


41/64 


145/128 




(0.5664) 


(0. 3301) 


(0.6406) 


(1.1328) 


3 


515/1024 


839/2048 


147/256 


687/512 




(0.5029) 


(0.4097) 


(0.5742) 


(1.3418) 


4 


0.4608 


0.4754 


0.5279 


1.5522 


5 


0.4298 


0.5328 


0.4927 


1.6834 



TABLE I: Hartree-Fock Quasiparticle (Aqp) and 
Skyrmion/Antiskrymion (Ask) particle-hole excitation 
gaps for the Dirac (D) bands of graphene and for parabolic 
bands (P). The two cases are distinquished by different 
dependences of form factor F n on Landau level (LL) index 
n. These results are for Coulomb interactions in an ideal 
two-dimensional electron system without finite thickness 
corrections and energies are in units of e 2 /el \Jn/2. In 
graphene the effective value of e depends on the dielectric 
screening environment provided by the substrate but is 
typically less than 2, whereas in GaAs and other common 
heterojunction systems e ~ 10. 



To determine whether or not skyrmions are the lowest 
energy charged excitation for a paraticular n we need to 
compare the energy of a skyrmion/antiskyrmion pair, 



8irp s 



(15) 



with the energy of an ordinary Hartree-Fock theory 
particle-hole pair, 



A PH = E{k -> oo). 



(16) 



Table [I] compares results for graphene Dirac-band and 
the ordinary parabolic band cases. In the Dirac-band 
cases both positive and negative values of n occur 
whereas the Landau level indices of parabolic systems 
are non-negative integers. In both cases, the increase in 
cyclotron orbit radius with |n| is reflected in the quan- 
tum form factor F n ; for the parabolic case the form fac- 
tor F n — L 2 l (q 2 £ 2 /2). Since the cyclotron orbit radius 
R c ~ (\/n it is clear simply on dimensional grounds that 
for Coulomb interactions and large n, Apn ~ e 2 /R c ~ 
e 2 /i x whereas Ask ~ e 2 R c /i 2 ~ e 2 /£ x ^/n. This 

difference in the large n behavior is already apparent in 
both cases in Table Q] For sufficiently large n then, 
Ask w ill exceed Aqp, the lowest energy charged exci- 
tations will be ordinary Hartree-Fock quasiparticles, and 
the ground state near integer filling factors will not have 
low-energy Skyrmion degrees of freedom. The quantita- 
tive calculations summarized in Table [I] show that or- 
dinary quasiparticles are already energetically preferred 
for ri = 1 in the parabolic band result obtained 

first by Wu and Sondhi, 1 - Interestingly the crossover to 
ordinary quasiparticles does not occur until n = 4 in the 
Dirac band case; thus skyrmion physics will occur within 
n = 0,n = l,n = 2 and n = 3 Landau levels in graphene. 
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We note in passing that we have so far compared only 
quasiparticle/quasihole and skyrmion/antiskyrmion pair 
excitation energies. In order to conclude that skyrmions 
and antiskyrmions are present in the ground state on 
both sides of integer filling factors, we need to demon- 
strate that the skyrmion energy and the antiskyrmion 
energy are separately smaller than the quasiparticle and 
quasihole energies respectively, whenever the pair exci- 
tation energies are smaller. In the case of N = 2 and 
v n = M = 1 this property follows^ from particle-hole 
symmetry. We demonstrate below that the pair excita- 
tion energy criterion also applies for graphcne, although 
the justification is subtly different for v n = 1 and v n = 3 
cases compared to the v n = 2 case. 

In earlier work Ezawa^fi and co-workers concluded that 
for v n = M > 1, lowest Landau level (LLL) skyrmions 
must have a charge that is a multiple of M. To address 
this discrepancy, we explicitly construct LLL skyrmion 
wave functions and demonstrate that it is indeed pos- 
sible to have charge ±1. Our microscopic single Slater 
determinant skyrmion wave function is constructed in a 
manner similar to that of Ref. 9. For dcfiniteness we 
discuss a hole-like skyrmion with q = — 1: 



AT* -2 



l^sky) = ["™ C M,m+l 



+ v mC M+lm 



m=0 

n 



n 



Ka<M 



|0>, 



(17) 



where the internal state labels match those of the symme- 
try broken ground state. Hartree-Fock skyrmion states 
are obtained by minimizing the expectation value of the 
Hamiltonian with respect to u m and v m . Skyrmion states 
with larger topological charges may be constructed in 
a similar manner. These single Slater determinant (or 
Hartree-Fock) skyrmion states correspond to classical 
skyrmions of the field theory l|12H , although the Coulomb 
self-interaction energy of skyrmions must be included in 
the field theory 8 in order to describe small skyrmions. 

The Hartree-Fock (or semi-classical) skyrmion state 
(I17|l has a rich internal structure, which may be ana- 
lyzed in a way similar to the analysis performed earlier 
on the ground state (JSJ. The state (|17|l is invariant only 
under unitary transformations among the first M — 1 
labels or the final N — M — 1 labels. It follows that 
the family of SU(N) transformations that actually trans- 
forms the skyrmion states form a coset space of dimension 
N 2 -(N-M-1) 2 -(M-1) 2 = 2M(N-M)+2(N-1). 
The set of transformations of dimension 2M(N — M) cor- 
responds to the order parameter fields of the ground state 
itself, and the additional 2(N — l)-dimensional space to 
(N — 1) skyrmion internal complex degrees-of-freedom. 
These must be specified in addition to location and size 
(determined by minimizing the energy with respect to 
u m and v m ) in order to completely characterize a classi- 
cal skyrmion. In other words, the skyrmion has N — 1 
internal (complex) degrees of freedom. A more intuitive 



way to understand this point is the following. Let us 
assume the ground state order parameter configuration 
is fixed (either spontaneously or by symmetry breaking 
pcrtubations). We now introduce a single skyrmion at 
the origin, and minimize its energy. This fixes the size 
of the skyrmion, and the magnitudes of u m and v m . On 
the other hand u m is a [CP(M — 1)] vector that lives in 
an SU(M) space spanned by the M occupied levels of 
the ground state, so there are M — 1 remaining internal 
(complex) degrees of freedom associated with it. Simi- 
larly there are N — M — 1 remaining internal degrees of 
freedom associated with v m . Adding the relative overall 
phase between u m and v m , we find the total internal de- 
grees of freedom is AT — 1. As a matter of fact u m and v m 
may be combined and viewed as an CP(N — 1) vector in 
the original SU (N) space; thus the N — 1 internal (com- 
plex) degrees of freedom of a skyrmion may be viewed as 
those of a CP(N — 1) superspini^i 

Once we move away from v n — M , we have a finite 
density of skyrmions in the ground state, provided that 
they are indeed the lowest energy charged excitations. 
For graphene we are able to judge the energetic ordering 
of charged excitations based on the energetic ordering 
of particle-hole excitations. To justify this statement at 
v n = M = 1, we must appeal to Zeeman coupling which 
is always present experimentally. Zeeman coupling se- 
lects a fully spin-polarized state and also selects fully 
spin-polarized skyrmions. We can therefore simply ig- 
nore the spin-degree of freedom and appeal to the same 
argumentaSi used for v n = M = 1 when N = 2. For 
v n = M = 3 we can then appeal to particle-hole symme- 
try in the N = 4 Landau level multiplet, which suggests 
that the situation is the same as v n = M = 1 . Finally for 
v n = M = 2 we can generalize the argument of Ref. 00] 
by appealing directly to particle-hole symmetry at this 
filling. 

In their classical ground state skyrmions will form a 
lattice. (Quantum corrections to the classical ground 
state become more important as the skyrmion density 
increases ii2i2S) In this case the internal degrees of free- 
dom associated with skyrmions on different sites will in- 
teract and the classical energy will be minimized by an 
arrangement with long-range order. When quantized, 
fluctuations in the internal degrees of freedom will re- 
sult in N — 1 spin wave-like modes, which will be present 
in addition to the single phonon mode associated with 
fluctuation in skyrmion positions. (In a strong magnetic 
field transverse and longitudinal position fluctuations are 
canonically conjugate leading 2 - to phonons with fc 3 / 2 dis- 
persion.) In the presence of full SU{N) symmetry these 
internal modes are gapless, and some (or all) of them may 
remain gapless in the presence of symmetry-breaking per- 
turbations, under appropriate conditions (see below for 
examples). For the SU(2) case we find N — 1 = 1 inter- 
nal mode, which is known previously^ here we provide 
a more general understanding of this result. 

For graphene N = 4, and the case v n = M = 2 is par- 
ticularly interesting. Weak Zeeman coupling will select 
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a unique fully spin-polarized, valley-singlet ground state. 
All collective modes are therefore gapped for u n exactly 
equal to 2, i.e. in the absence of skyrmions. In this case 
the SU (2) valley symmetry is preserved by the ground 
state. We nevertheless predict that the skyrmion lattice 
state near v n — 2 will have phonons and 3 additional 
gapless modes. Of these three internal modes, we predict 
that one has quadratic dispersion, while the other two will 
have linear dispersions. The quadratic mode may be un- 
derstood in the following manner. The skyrmion lattice 
state has a finite pseudospin (valley) magnetization, and 
is thus a spontaneous valley ferromagnet; the quadratic 
mode is simply the ferromagnetic spin wave in the valley 
channel. The two linear modes are the Goldstone modes 
corresponding to the two additional spontaneously bro- 
ken U(l) symmetries; these symmetries are broken by the 
fixed relative phases between u m and v m with the same 
valley index. Adding the phonon mode, we thus find one 
quadratic mode, one fc 3 / 2 mode, and two linear modes, 
all gapless. 

The situation is very different at v n = M = 1. In 
this case, the ground state is still fully spin-polarized 
due to the Zeeman coupling, but is also a spontaneous 
valley ferromagnet that breaks the SU(2) valley symme- 
try spontaneously. Thus we expect a single gapless mode 
with quadratic dispersion in the absence of skyrmions, 
the SU(2) valley pseudospin waves. The lowest energy 
skyrmion states are thus pseudospin textures in the valley 
degree of freedom^ and the skyrmion lattice states at 
low temperature will be fully spin-polarized. As a result, 
we expect only one gapless internal mode with quadratic 
dispersion which does not involve spin flips, while the 
other two internal modes which do involve spin flip will 
have a gap of gfisB. Thus the only new gapless mode 
of the skyrmion lattice state is the phonon mode in this 
case. 

Isolated individual quantum skyrmion states have well- 
defined total orbital angular momentum and internal 
SU(N) quantum numbers, which can be analyzed for 
v n = M = 1 following the strategy of Ref. [2CJ. These 
are the quantum counterparts of classical CP(N — 1) 
skyrmions. In this case the single Slater determinant 
state (0 is the exact ground state for 5-function interac- 
tion between electrons, since it has exactly zero energy. 
Hole-like skyrmion states may be identified as all zero en- 
ergy states for the case N$ = N e + 1. As emphasized in 
Ref. |2(] (see also Ref. l26j) . all such states may be written 
in the form 



(18) 



where the antisymmetric Jastrow factor Yii<j( z i~ z j) en ~ 
sures zero energy, while ^> b(z, x) is a bosonic wave func- 
tion that is symmetric under the exchange of spatial and 
internal coordinates of two particles. We can classify the 
skyrmion states based on the properties of this bosonic 
wave function. 



Due to the fact N$ = N e + 1, the bosons can only be in 
m = or m = 1 orbital states. Letting these occupation 
numbers be no and n\ respectively, the total angular mo- 
mentum of a state l|18[) measured from that of the ground 
state is AL = n%, thus no = N e ~AL. We now classify all 
the skyrmion states for a fixed AL based on the SU(N) 
representations they form (in the familiar SU(2) case 
these representations are labeled by a single quantum 
number, the total spin2£). Because bosons occupying the 
same orbital have totally symmetric orbital wave func- 
tions, their internal wave function must also be totally 
symmetric. It follows that the internal states of the no 
bosons in orbital m = form a totally symmetric repre- 
sentation of SU(N), as do the n\ bosons in orbital m = 1. 
In terms of Young tableau, they form two horizontal row 
representations, with no and n\ horizontal boxes respec- 
tively. More generally, each irreducible representation of 
the SU(N) group as represented by the Young tableau 27 
can be labeled by a set of N — 1 non-negative integers 
in descending order: [Zi , ^2 , • • • ,ijv-i], where each inte- 
ger represent the number of boxes in each row. Thus 
the two representations formed by the bosons in m = 
and m = 1 orbitals form representations [no, 0, • • • ] and 
[ni, 0, • • •] respectively. Now we take the direct product 
of these two representations and decompose them into 
irreducible representations of SU(N): 

[ n , 0, • • •] ® [m, Q, • • •] = [n + ni, 0, 0, • • •] 

© [n + ni - 1, 1,0, •••]©••• © [n o ,ni,0, • • •] 

= [N e , 0, 0, • • • ] © [N e — 1, 1, 0, •••]©•• • 

© [N e — AL, AL, 0, • • • ]. (19) 

In the above we have assumed that no > n±, or N e > 
2AL. If the opposite is true the positions of no and n\ 
need to be interchanged. The dimensionality of these 
representations may be found in the literature^ 7 .. This 
procedure classifies the hole-like skyrmion states at v n — 
1 based on their angular momentum quantum number 
and the irreducible SU(N) representations they form. 



IV. SU(N) FERROMAGNETS AT FRACTIONAL 
FILLING FACTORS 

We now turn our discussion to possible fractional quan- 
tum Hall (FQH) states (which have not yet been observed 
in graphene), many of which are also SU(N) ferromag- 
nets. Many FQH states may be constructed by starting 
with integer quantum Hall states^ and using the com- 
posite fermion (CF) flux-attachment ansatz. In this con- 
struction, we start from an IQH state with filling factor 
v n ,CF = fn, and attach an even 2n' flux quanta to the 
CFs to turn them into electrons; after the flux is spread 
out the electron filling factor in the valence LL becomes 

m 

= , T -i • (20) 
Zn'm ± 1 

This is, of course, the familiar Jain's sequence £££1 The 
difference here is that in the presence of the internal de- 
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generacy and SU (N) symmetry, the CF IQH state is an 
SU(N) ferromagnet as long as m is not a multiple of N. 
We thus expect the SU(N) symmetry properties to be 
reflected in the FQH states. The number of collective 
modes and the number of branches of skyrmion excita- 
tions should be the same as that of the IQH states at 
v n = M, if we identify M = mod(m, N). On the other 
hand if m is a multiple of N, we obtain an SU(N) sin- 
glet, and there is no spontaneous symmetry breaking; the 
system will be fully gapped. For cases with m < N, we 
write down below Laughlin-Halpcrin type of trial wave 
functions for the FQH states with the expected SU(N) 
symmetry properties in first quantization: 



ipn',m(z) 



N e 

m 

i<j 



2n 



a n n^- 

, a=lk<l 



Zl a ) 



Here A represents antisymmetrization of all coordinates, 
and we have neglected the common exponential factors 
for LL wave functions. The second factor is the first 
quantized wave function for J5|, while the first factor 
reflects flux attachment. The wave functions for m > N 
is more complicated, as in this case some of the CFs 
occupy higher LLs, and LLL projection is necessary^. 
We note that this type of construction was found to be 
generally reliable in predicting the spin structure in the 
SU(2) case without Zeeman splitting^. 

The low-energy physics of the FQH SU(N) ferromag- 
nets are also described by field theory (|12|l . In this case 
we do not have exact knowledge about the parameters 
(like the stiffness p s ) of the theory, but p 8 has been cal- 
culated numerically for the members of the Laughlin se- 
quence v = 1/3,1/5 for 1/r interaction based on the 
Laughlin trial wave function in the SU(2) casei^; they 
are 9.23 xKr 4 (e 2 /e£) and 2.34x l(T 4 (e 2 /ef) respectively. 
Using the same arguments as for the IQH case discussed 
above, we expect p s to take the same two values at the 
same filling factors in graphene. These values can be 
used to determine the energies of skyrmion-antiskyrmion 
pairs, which may be compared with transport measure- 
ments should FQH states be observed in graphene in fu- 
ture experiments. We caution however that at fractional 
filling factors, we can no longer appeal to particle-hole 
symmetry properties so that it is possible^ in general 
for Laughlin-like fractionally charged quasiparticles to be 
present in the ground state on one side of an incompress- 
ible filling factor and fractionally charged skyrmions on 
the other. 



V. CONCLUDING REMARKS 

In closing we comment on the possible effect of single- 
particle valley splitting on our St/ (4)-based analysis of 
graphene. We first note that the valley degeneracy of 
graphene is rather robust; for example unlike in Sili- 
con, simple strain cannot lift the degeneracy. Electron- 



electron interactions are likely to provide the most impor- 
tant source of Hamiltonian matrix elements that break 
the valley portion of the ST/ (4) symmetry. Our consider- 
ations should nevertheless largely apply for weak valley- 
symmetry breaking, as long as the characteristic energy 
scale of these terms is much weaker than the Coulomb 
interaction scale. We mention that among the possible 
extrinsic sources of 4-fold degeneracy lifting in the single- 
particle Hamiltonian are edge effects and inter-valley 
scattering. For large extrinsic splittings, the SC/(4) sym- 
metry will be reduced to SU(2) xU(l)x U(l), where the 
remaining SU (2) symmetry corresponds to spin, which is 
further reduced by Zeeman splitting as discussed earlier. 
At the lowest temperatures and energy scales the physics 
may be be similar to that of a bilayer system^, which 
may support other types of broken symmetry state a 32,33 . 
The question of whether our predicted interaction-driven 
spontaneous breaking of valley degeneracy is playing the 
key role in the observed valley splitting at high field 4 
can be decided by careful measurements of zero-field val- 
ley splitting. If the experimental zero-field valley split- 
ting is negligibly small, then it seems certain that SU (4) 
quantum Hall ferromagnetism, associated with the spon- 
taneous breaking of valley and spin degeneracy, is al- 
ready playing a role in the high-field quantum Hall ex- 
periments in graphene. (It seems clear that many-body 
physics does not have a large influence on quasiparticlc 
valley-splitting in the absence of a field.) Direct observa- 
tions of skyrmions and associated collective excitations in 
graphene then take on particular experimental relevance. 
We further note that both Zeeman and valley splittings 
are single electron effects; electron-electron interactions 
also have weak symmetry breaking effects due physics at 
lattice scale, as discussed recently^ 

In summary, we have identified exact broken sym- 
metry eigenstates of S[/(A r )-invariant Hamiltonians in 
the quantum Hall regime of graphene. We argue that 
these states are ground states for any physically sensi- 
ble repulsive interaction between the particles, and for 
Coulomb interactions between electrons in particular. 
Given SU(N) invariance we were able to obtain exact 
results for the elementary collective excitation spectrum. 
We found that in the absence of symmetry breaking there 
are M(N ~ M) gapless collective modes when M of the 
TV members of a Landau-level multiplet are occupied. In 
addition, we have shown that for Coulombic electron- 
electron interactions, skyrmions are much more robust in 
systems like graphene with a Dirac band structure, than 
for the more familiar parabolic band systems. We pre- 
dict that dense-skyrmion systems will be ubiquitous in 
the quantum Hall regime of graphene and that skyrmion 
lattice states near even integer filling factors will have 
four branches of gapless collective modes, two with k l 
dispersion, one with fc 3 / 2 dispersion and one with k 2 
dispersion. The identification of probes that can study 
skyrmion physics in graphene layers is therefore an at- 
tractive challenge for experiment. 
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